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Abstract 

. We continue our investigation of the nested loop approach to the 0(n) model 

on random maps, by extending it to the case where loops may visit faces of 
arbitrary degree. This allows to express the partition function of the 0(n) loop 
model as a specialization of the multivariate generating function of maps with 
controlled face degrees, where the face weights are determined by a fixed point 
condition. We deduce a functional equation for the resolvent of the model, 
involving some ring generating function describing the immediate vicinity of 
^ ; the loops. When the ring generating function has a single pole, the model is 

amenable to a full solution. Physically, such situation is realized upon consid- 
ering loops visiting triangles only and further weighting these loops by some 
local bending energy. Our model interpolates between the two previously 
solved cases of triangulations without bending energy and quadrangulations 
with rigid loops. We analyze the phase diagram of our model in details and 
derive in particular the location of its non-generic critical points, which are 
in the universality classes of the dense and dilute 0(n) model coupled to 2D 
quantum gravity. Similar techniques are also used to solve a twisting loop 
model on quadrangulations where loops are forced to make turns within each 
visited square. Along the way, we revisit the problem of maps with controlled, 
possibly unbounded, face degrees and give combinatorial derivations of the 
one-cut lemma and of the functional equation for the resolvent. 
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1. Introduction 



The study of planar maps, i.e. proper embeddings of graphs in the two-dimensional 
sphere, has been the subject of an intense activity in combinatorics since the seminal 
papers of Tutte in the 60's Q. Recent developments @ entered the realm of probability 
theory, by the study of particular ensembles of maps drawn according to some prescribed 
probability distribution. Such ensembles arise naturally in physics where maps serve 
as discrete models for various types of fluctuating surfaces such as fluid membranes in 
soft condensed matter physics or quantum space-times in the theory of two-dimensional 
quantum gravity. In these contexts, maps have often been equipped with additional 
degrees of freedom such as spins or particles, giving rise to nice universal critical be- 
haviours in the limit of large maps || . A particular important class of such models are 
the so-called 0(n) loop models, where maps carry self- and mutually-avoiding loops, 
each weighted by n [|]-|10]]. New universality classes are reached in these 0{n) loop 
models when the loops become large and modify the statistics of the underlying map. 
While the critical exponents have been determined exactly and corroborate the cele- 



brated KPZ relations which relate them to critical exponents of the 0(n) model 
on a fixed regular lattice, little is known on the random geometry of these models. 

This has to be contrasted with the case of planar maps without additional degrees 
of freedom, whose geometry is now well understood. In particular, when maps are 
endowed with the graph distance, thus are viewed as metric spaces, their generic scaling 
limit (in the Gromov-Hausdorff sense) is described by a unique remarkable probabilistic 
object: the Brownian map |2,12,13|. Physically, it corresponds to the universality class 



of so-called pure gravity. One may escape from this universality class by considering 
maps with large faces (or "holes") ||14|| . More precisely, one works with the so-called 
Boltzmann ensemble of random maps, where each face receives a weight depending 
on its degree. If the face degrees are bounded (i.e. the weight for faces of degree k 
is set to zero for k large enough - which is a usual assumption), the scaling limit is 
the Brownian map. Non-generic scaling limits may only be obtained by allowing for 
unbounded degrees, and furthermore fine-tuning the weights in such a way that the 
degree distribution of a typical face has a heavy tail characterized by an exponent a G 
]1,2[. Under this assumption, the scaling limit is a so-called "stable map" of Hausdorff 
dimension 2a G]2,4[ [T3], different from the Brownian map of dimension 4. At first, 
the fine-tuning of the weights might seem slightly unnatural, but it was proposed that 
it occurs spontaneously in the context of critical 0(n) loop models on random maps: 
considering a sample configuration of such a model, if one erases all the outermost loops 
and their contents, the resulting map (called gasket) is a map drawn according to a 
Boltzmann ensemble with non-generic scaling limit. 

This mechanism was investigated in a recent paper |15|], for a particular class of 
0(n) loop models on tetravalent maps (dual to quadrangulations) . Notably, n was 
shown to be generically related to a via n = 2 sin ira (where n should vary between 
and 2 in order for the loop model to have critical points, the two corresponding possible 
values of a being associated respectively with the dense and dilute universality classes) . 
A full explicit solution of the model was obtained for a specific model, the so-called 
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rigid loop model where loops are forced to go straight within each visited square. The 
purpose of the present paper is to extend the study of ]15| to the largest possible class 
of 0(n) loop models amenable to a full solution by the same techniques. As we shall 
discover, this class corresponds to models with loops visiting triangles only, and with 
some additional bending energy term which controls the rigidity of the loops. These 
models interpolate, upon varying the bending energy, between the standard 0(n) loop 



model on triangulations j§-[10| an d the rigid loop model of |L5| on quadrangulations 
(with squares formed of two triangles). Note that here, we define our model only on 
planar maps. Extending and solving it for maps of arbitrary topology would present 



no difficulty, following the lines of |1J|, although we do not address it here. Indeed, 
once the planar case is solved, the answer for other topologies should be given by a 
topological recursion formula. 

The paper is organized as follows: in Section 2, we recall the nested loop approach 



initiated in |T5| as a general strategy that allows to reformulate any 0(n) loop model 
as a model of maps with controlled face degrees. In practice, the nested loop approach 
consists of a gasket decomposition (Section 2.1) obtained by cutting the maps along 
the superimposed loops, resulting in a bijective coding of the 0(n) loop configurations. 
This coding translates into a fixed point condition for the degree-dependent weights in 
the equivalent problem of maps with controlled face degrees (Section 2.2). This in turn 
translates into a functional equation for the so-called resolvent of the 0(n) loop model, 
which is the generating function for loop configurations on maps with a boundary of 



controlled length (Section 2.3). Here, in contrast with \TB\, we use a shortcut to directly 
write down the functional equation from the fixed point condition. Explicit examples 
of functional equations are given in Section 2.4 for a number specific models that have 
been considered before. In all generality, the functional equation involves a bivariate 
ring generating function which is the (grand canonical) generating function for the im- 
mediate surroundings of the loops. Section 3 is devoted to finding the largest possible 
class of models for which the functional equation may be solved by a straightforward 



generalization of the techniques presented in |T5| and originally developed by Rostov 



Rostov and Staudacher 0, Eynard and Zinn-Justin |7j and Eynard and Rristjansen 
§, and reformulated later in terms of algebraic geometry ||16|1 . These models corre- 



spond to having a bivariate ring generating function with a single pole, as explained 
in Section 3.1, in which case the functional equation simplifies drastically and involves 
crucially some (general) decreasing homographic involution. Interestingly enough, all 
these models may be realized by considering a model of loops visiting triangles only and 
weighted by some particular bending energy, as explained in Section 3.2. The details 
of the solving technique are finally recalled in Section 3.3. This strategy is applied ex- 
plicitly in Section 4 to study the phase diagram of the 0{n) loop model with bending 
energy defined on triangulations, i.e. when the unvisited faces are themselves triangles. 
An explicit derivation of the non-generic critical line, corresponding to the location in 
the parameters of the model of its dense and dilute phases, is given in Section 4.1 while 
the whole phase diagram is discussed in Section 4.2. The limit n — > of the model is 
discussed in Section 4.3 as a number a combinatorial simplifications arise in this case. 
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Section 5 is devoted to the solution by similar techniques of the so-called twisting loop 
model on quadrangulations which is yet another particular case of the model considered 
in [13] where, in contrast with the rigid loop model, loops are forced to make a turn 
inside each visited square. Our derivation of the functional equation for the resolvent of 
the 0(n) model relies on the analytic properties of the similar resolvent in a model of 
maps with controlled face degrees. We review these properties in Section 6. We present 
in Section 6.1 a new combinatorial proof of the so-called one-cut lemma and rederive in 
Section 6.2 the functional equation satisfied by the resolvent. Section 6.3 is devoted to 
the proof of a technical result while the case of complex-valued face weights is briefly 
discussed in Section 6.4. We gather our conclusions in Section 7. 



2. The nested loop approach 

2.1. The gasket decomposition 

In this paper, we consider loop configurations on planar maps defined as follows: 
given a map, a loop is an undirected simple closed path on the dual map (i.e. it visits 
edges and vertices of the dual map, hence visits faces and crosses edges of the original 
map). A loop configuration is a set of disjoint loops. Choosing a loop configuration 
amounts to choosing a set of crossed edges in such a way that each face is incident to 
exactly or 2 crossed edges (with multiplicities, to account for self- folded faces). We 
call map with a boundary of length p {jp > 1) a rooted planar map where the external 
face (i.e. the face on the right of the root edge) has degree p. By convention, when 
considering loop configurations on maps with a boundary, we demand that the external 
face is not visited by a loop. The external face allows to distinguish the exterior and 
interior of a loop. In particular, we may associate to each loop its outer (resp. inner) 
contour consisting of all the exterior (resp. interior) edges that are incident to a face 
visited by the loop. Each contour is a closed, but non-necessarily simple, path on the 
map. 

Following [|I5J , let us define the gasket decomposition of a map with a boundary of 
length p endowed with a loop configuration. It consists in cutting the map along the 
outer and inner contours of each outermost loop. This splits the map into connected 
components which are: 

- its gasket which is the map spanned by the edges which were exterior to all the 
loops, 

- for each outermost loop: 

- its ring formed by the faces visited by the loop, 

- its internal map which is the map spanned by the edges which were interior to 
the loop. This map is endowed with a loop configuration consisting of all the 
loops which were originally interior to the outermost loop at hand. 

All these components are rooted canonically, as displayed in Fig. 1. Thus, each internal 
map may be viewed as a map with a boundary of the same nature as the original map. 

By construction, the gasket is itself a map with a boundary of length p (note that 
the gasket contains all the edges incident to the external face), without loops, and whose 
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Fig. 1: A schematic picture of the gasket decomposition. Cutting the map 
along the outer contour of each outermost loop and removing the interior of 
these outer contours yields the gasket, with a new type of faces: the holes. 
Upon cutting along the inner contours of the same loops, the former content 
of each hole may itself be divided into a ring and an internal map. All these 
components are rooted by some canonical (yet irrelevant) procedure. We 
may for instance consider the leftmost loop-avoiding shortest path from 
the root edge of the map (i.e. the root of the external face) to the outer 
contour of any outermost loop. This selects an edge on this contour which 
can be used as root for the corresponding ring. By then transferring the 
root from the outer to the inner contour by some local rule, it gives a root 
for the corresponding internal map. 

inner faces are of two types: regular faces corresponding to faces of the original map 
that were exterior to all the loops, and holes delimited by the outer contours of the 
former outermost loops. Note that an outermost loop having outer and inner contours 
of respective lengths k and k' gives rise to a hole of degree k in the gasket and to an 
internal map with a boundary of length k' . Furthermore, both the gasket and internal 
map may contain separating vertices. In particular, a multiple point along the outer 
(resp. inner) contour yields a separating vertex incident to a hole in the gasket (resp. 
to the external face of the internal map). By contrast, we decide that multiple points 
on the contours are split into distinct vertices in the ring, so that both sides of the ring 
are simple paths of lengths k and k' . 

To summarize, the gasket decomposition provides a bijection between, on the one 
hand, maps with a boundary of length p endowed with a loop configuration, and on 
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the other hand, maps with a boundary of length p with regular faces and holes, where 
each hole of degree k (k > 1) is endowed with (i) a ring with sides of lengths k and 
k' for some k! > and (ii) a map with a boundary of length k' endowed with a loop 
configuration (by convention there is a unique map with a boundary of length 0, the 
vertex-map, with no loops). 

2.2. The fixed-point condition 

We may use the above gasket decomposition to relate the partition function of 
various loop models on random maps to the multivariate generating function for maps 
with prescribed face degrees. More precisely, consider an 0(n) loop model on maps 
with a boundary where the Boltzmann weight of a configuration (i.e. a map endowed 
with a loop configuration) is the product of: 

- for each k > 1, a weight per inner face of degree k not visited by a loop, 

- a weight per loop that only depends on the structure of its ring, that we write in 
the form nw r where n is the loop fugacity and w r is, at this stage, an arbitrary 
function of the ring r. In all the cases discussed in this paper, w r will be given as a 
product of some local weights, leaving n as the only non-local weight in the model. 

Let us denote by F p the partition function for such an 0(n) loop model on maps with a 

boundary of length p, which depends on n, the jf''s and the w r 's (by convention, we set 
Fq = 1). It follows from the gasket decomposition that F p is obtained by specializing a 
multivariate generating function for maps without loops but with controlled face degrees. 
More precisely, we denote by J- p (gx, g%, . . .) the multivariate generating function for maps 
with a boundary of length p with weight g k per inner face of degree k (by convention 
J-q = 1), and we introduce the ring generating function 

A k ,k' = ^2 Wr ( 2-1 ) 

rETZ(k,k') 

where TZ(k, k') is the set of rings with outer length k and inner length k', rooted on the 
outer contour. Then, the gasket decomposition translates into the relation 

F p =T p ( 9l ,g 2 ,...) (2.2) 
at values of the weights g k satisfying the fixed-point condition 

9k = 9 k + n A k>k/ Tk' (ffi, t/2, • • •) for all k > 1 . (2.3) 

fc'>0 

Let us now give a few examples of ring generating functions. If we assume that 
the loops only visit triangles, and that the weight for a ring r made of m triangles is 
w r = h m , then we recover the 0{n) loop model studied in |5f||. We have 

<V = (" + ^ _1 )^ +fc/ , (2-4) 



6 



as seen by the following simple counting argument: a ring with outer length k and inner 
length k' is a sequence of k + k' triangles, k of which are facing outwards, that we may 
read starting from the distinguished outward-facing triangle carrying the root of the 
ring. 

Consider now the 0(n) loop model studied in [15| . Here, the loops only visit 
squares, and there are two types of visited squares: those such that the loop crosses 
opposite sides (type b) and those such that the loop crosses adjacent sides (type c). The 
weight for a ring r made of m squares of type b and m' squares of type c is w r = h™h™ , 



and the ring generating function is given by flT5| , eq. (2.4)] 



min(fc,fc') 
a quad _ ^ 
A k,k' ~~ 



2k 



(k + k')\ 



3=0 



k + k> (2j)\(k-j)\(k>-j)\ 



hl j h 



k+k'- 

2 



2 J 



(2.5) 



In particular, for h 2 = we recover the so-called rigid loop model for which 

Kf = ^h\ (2.6) 

as seen directly by noting that the unique ring of outer length k is made of k squares 
and has inner length k! = k. Alternatively, let us consider the opposite case h\ = that 
we call the twisting loop model. Clearly, the associated ring generating function A^^ 
vanishes whenever k or k' is odd, and 

A^%,=^ + l~ l )h k 2 +k \ (2.7) 

as seen directly by a counting argument similar to that for (|2.4|) (there are 2 possible 
choices for the position of the root on its incident outward-facing square, thus the extra 
factor 2). 

In general, we may encode the A/^/'s into the grand generating function 



A( X} y) = j2J2 A ^ ,xk y k ' ■ ( 2 - 8 ) 

k>l k'>0 

The apparent asymmetry between k and k! in the above expressions is related to the 
fact that the ring has a distinguished edge on the outer contour. Assuming that the 
ring weight does not depend on the position of the distinguished edge, we may write 

d 

A(x,y) = x—\ogH(x,y) (2.9) 

where log H (x, y) is now the generating function for unrooted rings (involving symmetry 
factors). If the ring weight is invariant when exchanging the exterior and the interior, 
then H (x, y) is symmetric in x and y. In all the cases discussed in this paper, H (x, y) 
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thus A(x, y) will be simple rationals function of x and y. For instance, in the triangular, 
quadrangular, rigid and twisting cases discussed above, we have respectively 



1 — h(x + y) " 1 — h(x + y) 

H^ ad (x,y) = — 1 —— 9 A^ d (x,y) h lX y + 2h 2 x 2 



1 — h\xy — h2(x 2 + y 2 ) ' " 1 — h\xy — \i2{x 2 + y 2 ) 

H r ^ d (x, y) = , A^ id (x, y) = : klX h V , 

1 — h\xy 1 — h\xy 

H twiS \x, y) = : r A , or , A twist (x, y) 2/ ' 2 -''"' 



1 — /i2(a; 2 + y 2 ) ' ' 1 — h 2 (x 2 + y 2 ) 

(2-10) 

More generally, if the ring weight is simply given by a product of local face weights 
(Ifw), say hi^i per face incident to I external edges and I' internal edges, then H(x,y) 
corresponds to the generating function for unrooted rings with a distinguished face. 
Such rings are in correspondence with sequences of faces, so that 

which are rational if the family (/i^')^,£'>o has finite support. The symmetric case 
corresponds to having hi,i> = for all £, £' > 0. 

2.3. The functional equation for 0(n) loop models 

The fixed point condition ( [2.3| ), which is an equation for the infinite sequence 
(gk)k>i, ma y be rephrased as a functional equation for the resolvent of the model, 
defined as 

p>0 

Recall that F p is a specialization of the generating function J- p (gi, g 2 , ■ ■ •) for maps 
with controlled face degrees. As such, the resolvent W{x) satisfies a number of known 
analytic properties, which we review in Section 6. In particular, when the model is 
well-defined, that is to say the parameters of the model (n, the s and the uv's) are 
such that the F p 's are finite, W may be analytically continued into a function that is 
holomorphic on the complex plane except on a cut where it has a finite discontinuity 
(see Section 6.1). When all weights are real non-negative, this cut is a real interval 
[7_,7_i_] with 7_|_ > |7_|. The discontinuity of W on its cut is the so-called spectral 
density 

= w(x . i0):W(x + i0) 

2m 

which vanishes at x = y±. Moreover, as discussed in Section 6.2, we have the funda- 
mental relation 

W(x + i0) + W(x-i0) = V'(x) , xe]7_,7+[, (2.14) 



8 



where we introduce the shorthand notation 

V'(x)=x-J29kX k - 1 . (2.15) 



k>i 



This fundamental relation holds for the resolvent of general maps with controlled face 
degrees, i.e. does not require that the weight sequence (gk)k>i satisfies the fixed-point 
condition ( |2.3| ). This latter condition allows to rewrite V'(x) as 



k>l k>l k>>0 , n 1 -v 

- ( 2 - 16 ) 
= V£(x)-—<b A(x,y)W(y)dy 

where A(x,y) is defined as in (|2.8|), C is a contour surrounding [7_,7 + ] and 

V^(x) = x-J29 { k ° ) x k - 1 . (2.17) 

k>l 

Note that the contour C must be included in the domain of analyticity of y \-> A(x, y) 
in order for V'(x) to be well-defined. We end up with the linear functional equation 

W(x + iO) + W(x - iO) = V&x) " jf A (x, y)W(y)dy , x e] 7 -, 7+[ ■ (2-18) 

A new feature of this equation is that W now appears also on the rhs. 

2.4- Some examples of functional relations 

When A(x, y) is a rational function, we may evaluate the integral in ( |2.16| ) using 
the residue theorem. For instance, in the triangular case, A tri (x , y)W (y) has by ( |2.10|) 
one simple pole at y = h~ l — x, hence we recover the functional equation [|]-[| 

W tri (x + iO) + W tri (x - iO) = Vl{x) - n W tYi (h- 1 - x). (2.19) 

Similarly in the rigid case, A rigld (x, y)W(y) has one simple pole at y = l/(h\x) and a 



residue 1 at infinity (since W(y) ~ 1/y), hence we recover the functional equation [T5 



W rig[d (x + iO) + W lig[d (x-iO) = VMx) + - - -^r igid ( V-) • (2.20) 

x h\x z \h\x J 

In the generic quadrangular case Q12 > 0), A quad (x, y) W(y) has two simple poles at the 
roots 

y ± (x) = Z^Wg|Ejgg±gg (2 . 21) 
9 



and the functional equation reads 

W quad (x + iO) + W quad (x - iO) = 

(2 22) 

Vb'(x) +ny' + (x) W quad (y+(x)) +ny'_(x) W quad (y-(x)) . 

Let us now discuss the case where the gasket is a bipartite map, i.e. when all its 
faces have even degree. Such a situation occurs whenever g^k+i an d A 2 k+i,2k' vanish 
for all k, k' > 0, so that g 2 k+i = for all k. Then F p vanishes for all p odd so that W (x) 
is an odd function of x (in particular, 7_ = —7+). We may then write W(x) = x W(x 2 ) 
where VF(X) has a cut [0, T] (with r = 7? = j_) on which it satisfies 

W(X + i0) + W(X -iO) = V'{X), Xe[0,T[ (2.23) 



with V'(X) = 1 — J2k>i 92kX k 1 . This equation is the bipartite counterpart of ( |2.14| ). 
The fixed point condition (|2.3|) then yields the functional equation 



W(X + iO) + W(X - iO) = Vq{X) - I A{X, Y)W(Y)dY (2.24) 

2i7rA J 

where the integral is over a contour encircling the cut [0, V] and 

Vq{X) = 1 -Y,9 i 2 k ) X k - 1 , A(X,Y) = J2J2 Mk,2 k >X k Y k ' . (2.25) 

k>l k>lk'>0 



In particular, in the twisting loop model, A tw[st (X, Y) = A tw[st (VX, y/Y) has by ( gig ) 
a simple pole at Y = h^ 1 — X with residue 2 so that the functional equation reads 

W(X + iO) + W(X - iO) = Vq(X) - 2n W(h^ - X) , (2.26) 



which is very similar to ( |2.19|) . The solution of the triangular case (|2.19|) was the topic 



of the original work of Rostov 0, and Eynard and Kristjansen [@,|9|]. With identical 
techniques, the solution for the twisting loop model on quadrangulations ( |2.26|) can be 
easily deduced (see Section 5). With a few modifications, the method was adapted to 
solve ( p.20| ) in ||15|| . In this article, we show how to solve via similar methods a more 
general family of models (which contains those three examples as special cases). 



3. The 0(n) loop model on cubic maps with bending energy: equations 

3.1. Ring generating functions with one pole 

A peculiar feature of the triangular and rigid loop models is that both the ring 
generating functions A tri (x,y) and A Tlgld (x, y), as defined in ( |2.10| ), are rational func- 
tions with one pole, which leads to a simplification of the functional equation (j2.18|) 
into respectively ( |2.19| ) and ( |2.20|) . In this section, we put those two models under 
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the same roof by considering the general "one-pole case" (note that the bipartite ring 
generating function A twist (X, Y) is also a one-pole rational function but, for simplicity, 
we postpone the discussion of the twisting loop model to Section 5). 

More precisely, let us assume that the ring weight is symmetric and does not depend 
on the position of the distinguished edge, so that A(x, y) is given by (|2.9|) , with H(x, y) 
a symmetric one-pole rational function of the form 

H{x,y) = ^ r (3.1) 
x - s(y) 

Note that the numerator cannot depend on x, since any zero would give rise to a pole 
in A(x, y). Demanding that H(x,y) be symmetric in x,y yields a functional equation 
for h and s. It is readily solved (e.g. write the equation for y = and solve for h(x), 
then substitute into the equation for y = 1) and implies that s(x) takes the form 

•W " (3 ' 2) 

for some real parameters a, (3 and 5, while h(x) oc If {a — f3x). Note the remarkable 
property s(s(x)) = x. In other words, we recognize in ( |3.2| ) the general expression for a 
homographic involution. In particular, this involution reduces to s(x) = h~ x — x for the 
triangular model and to s(x) = l/(h\x) for the rigid model. By (|2.9|), we have 



A(x,y) = - r ^ (3.3) 

s{y) - x 

hence the ring generating function Ak,y satisfies the exponentiation property 

A k ,k>y k ' = s(y)-\ k>l. (3.4) 

k'>0 

For positive ring weights, s(y) is then necessarily a decreasing positive function of the 
real positive variable y. This imposes s(0) > together with the extra condition (3 2 — 
ad > 0. Hence, the involution has two real fixed points at (/3 ± a//3 2 — a5)/5. We will 
see in the next subsection that all the decreasing homographic involutions are indeed 
obtained in a simple model of loops with bending rigidity. In Appendix A, we mention 
the cases where the homographic involution is increasing, although they cannot be 
reached with non-negative local weights. 

By the form ( |3.2|) for s(x), one may easily check that A(x,y), as given by (|3.3|) , 
may be written in the alternative form 

M*,V) = J ?^ + ? £M, (3-5) 
y — s(x) 2s'{x) 

which displays its single pole in the variable y at y = s~ 1 (x) = s(x). We may now 



evaluate the contour integral in (|2.18| ) via the residue theorem. Recall that, in order 
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for the model to be well-defined, the cut [7_,7 + ] of W has to be contained, for all x 
on the cut, within the disk of convergence of y i— > A(x,y) which here has radius |s(x)|. 
This shows that the interior of the cut may not overlap with its image under s (but the 
endpoints of the cut may be fixed points of s). The integrand A(x, y)W(y) has a pole 
at y = s(x), with residue equal to —xs'(x)W(s(x)). Since W(y) ~ l/y for y — > oo, the 
residue at infinity is equal to lim y _ 5 . 00 A(x,y) = xs"(x)/(2s'(x)). We deduce that 

i j A(x, y)W(y)dy = -xs' (x)W(s(x)) + (3.6) 

which, inserted in (|2.18|) , yields the functional equation 

W (x + iO) + W (x - iO)- ns' (x)W (s(x)) = V^(x)- —- ^ , xe] 7 _,7 + [ (3.7) 

This equation generalizes eqs. ( |2.19|) and ( |2.20|) , which are recovered for s(x) = h~ x —x 
and s(x) = l/(h\x) respectively, as expected. 

3.2. Realization in a model of loops with bending energy 

In this Section, we exhibit a model whose ring generating function reproduces pre- 
cisely the most general symmetric one-pole rational function of the previous subsection, 
hence whose resolvent is the solution of the general functional equation (|3.7|) . We con- 
sider again an 0(n) loop model where the loops visit only triangles, with a weight h per 
visited triangle, but we also introduce a bending energy for the loops in the following 
way: each triangle visited by a loop has exactly one edge which is not crossed by the 
loop. We say that the triangle faces outwards (resp. inwards) if this edge is on the same 
side of the loop as the external face (resp. on the opposite side). Now for each pair of 
successive triangles along a loop, we attach a weight a if the triangles are of the same 
nature (either both facing outwards or both facing inwards) and a weight 1 otherwise. 
We may interpret — log a as a bending energy for the loops as this new weight favors 
loop turns for a > 1 and penalizes them for < a < 1. Such a bending energy was 
previously introduced in JT^] in the different context of Lorentzian triangulations. Note 



that our model precisely interpolates between the regular 0(n) loop model on triangles, 
recovered for a = 1, and the rigid 0(n) loop model on squares, recovered for a = and 
h± = h 2 . Indeed, at a = 0, outward- and inward- facing triangles necessarily alternate 
along the loop and, upon concatenating them by pairs, we may simply replace them by 
squares (with weight h?) having opposite sides crossed by the loop. 

Let us now compute the ring grand generating function A(x, y). Rings whose inner 
contour has zero length are non-empty sequences of outward-facing triangle only, each 
triangle receiving a weight ahx. These rings therefore have a grand generating function 
equal to tztjt- ■ Rings with inner contour of non-zero length may be viewed, before 
choosing their root edge, as cyclic sequences of composite objects consisting of (i) a 
non-empty sequence of outward-facing triangles followed by (ii) a non-empty sequence 
of inward-facing triangles. Each composite object therefore yields a grand generating 
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function 1 ^ hx ■ i^hy ( w ith no weight a for the first triangle of each sequence since it 
follows a triangle of the opposite type). We deduce eventually 

ahx d ( , / hx hy \ \ . 
A&y)=i r + ^ -log 1-- =- - \ (3.8) 



1 — ahx dx \ \ 1 — ahx 1 — ahy 

where, in the second term, the log(-) accounts for the cyclic sequence of composite 
objects while the xd x accounts for the rooting of the ring. We observe that A(x,y) is 
indeed of the form ( |2.9| ) with 

^ 1 — ah(x + y) — (1 — a 2 )h 2 xy ^ ^ 

which is a symmetric one-pole rational function (observe the similarity with Eq. (2.3) 
of fT^]). We are therefore precisely in the situation discussed in Section 3.1, hence we 
may rewrite A(x, y) in the forms (|3.3|) and (|3.5| ) with s(x) now given by 



/ \ 1 GbJ~b3C , v 

six) = — 9 , ; » (3.10) 

v ' ah + (l- a 2 )h 2 x v ' 



which has two fixed points at x = ( a ±^ h ■ We therefore recover the general form (|3.2|) for 



an arbitrary decreasing homographic involution such that s(0) > (with = ^ = ha, 

— = {a 2 — l)h 2 , (3 2 — a5 = h 2 a 2 > 0). In other words, our loop model with bending 
energy realizes precisely the general "one-pole" case. Note that s(x) reduces to h~ x — x 
for a = 1 and to -A- for a = 0, as it should. 

3.3. Solving strategy 

The technique developed by Eynard and Kristjansen || to solve the functional 
equation (|2.19| ) may be adapted to ( |3.7| ) with few modifications, as we will now sketch. 
Briefly, we may see that, for fixed given y±, the functional equation admits a unique 
solution W analytic outside the cut [7_,7 + ] and such that W(x) ~ 1/x for x — > oo. 
Then, the unknowns r )± are determined a posteriori by the condition p( / y±) = 0, where 
p is the spectral density ( |2.13|) . 

We will restrict ourselves to polynomial Vq, i.e. we assume that the degrees of the 
faces not visited by loops are also bounded. In that case, Vq is defined on the whole 
complex plane, and we easily find a particular solution of (577) with no discontinuity 
along the cut, namely 

_ 2Vi(x) + ns\x)V>(s(x)) ns"(x) 
WpaAx} ~ 4^ 2(n + 2)s'(xY (3 ' n) 

To check that this is indeed a particular solution, note that, s(x) being an involution, 
we have s'(s(x)) = l/s'(x) and s"(s(x)) = —s"(x)/s'(x) 3 . Since Vq is a polynomial, 
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observe that W par t has (multiple) poles at oo and at s(oo). Any solution of ( |3.7| ) is of 
the form 

W{x) = W part (x) + W hom (x) , (3.12) 
where Whom is now a solution of the homogeneous functional equation 

W hom (x + iO) + W hoia (x-iO)-ns'(x)W hom (s(x)) =0 , xG] 7 -,7+[. (3.13) 

Assuming that s(oo) 7^ 00, viz. 5 7^ in ( |3.2|) and a 7^ 1 in ( |3. 10|) , the requirement 
W(x) ~ 1/x for x — > 00 yields the condition 

W hom (x) = -W part (z) + I + O ( ±) = + + O (\). (3.14) 

x \x A J 4 — n z n + 2x \x A J 

Moreover, the analyticity of W outside the cut, or the finiteness of p on the cut, requires 
that the pole of W part at s(oo) is canceled, which may be rephrased as 



s'(x)W hom (s(x)) = -s\x)W p&rt (s(x)) + o(-^j 



nVfa) n 1/1 



4 — n 2 n + 2 x \x 2 

^W hom (x) + O ( 
2 \x z 



(3.15) 

for x 00. When s(oo) = 00, viz. 5 = in ( |3.2|) or a = 1 in ( |3.10| ), the two conditions 
( |3.14[ ) and ( |3.15|) are replaced by the single condition 



= + i + o (±) = + i + o (i 

(3-16) 

The general solution of ( |3.13| ) may be expressed in terms of elliptic functions: let 
us introduce the elliptic integral 

V[x) = I V±(£ - 7+)« - 7-)« - *(7+))(* ^(r0) ' (3 ' 17) 

Since s is (locally) decreasing and the cut [7- , 7+] does not overlap with its image, 
the sign may be chosen in such a way that the polynomial under the square root be 
positive over [s( 7 +), s( 7 _)] (when 5(7+) > s(7_), this interval is to be understood as 
the union ] — oo,s( 7 _)] U [s( 7 +), oo[). In particular, this implies that T = i?(s( 7+ )) is 
real and positive. By the same reason, the polynomial under the square root is negative 
on [s(7_),7_] (understood as ] — 00, 7_] U [s( 7 _),oo[ if 7_ < s( 7 _)) so that i>(7-) 
has two pure imaginary determinations ±iT" (T' > 0). Finally, we may check that 
v ( 7+ ) = T±iT'. 

For the moment, we assume that s(p(±) 7^ r f± so that T < 00. The function 
x 1 — y v(x) has branch points of order 2 at 7±,s(7±) and maps the lower (resp. upper) 
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half-plane onto the rectangle [0, T] +i[0, T'] (resp. [0, T] +i[— T f , 0]). Its inverse function 
v i — y x(v) may be analytically continued to an even doubly periodic function with periods 
2T and 2iT', with simple poles at ±t>oo mod 2TZ + 2iTZ, where = lim x ^ 00 v(x). 
By the change of variable x — > s(x) in (|3.17| ), we readily see that 

x(v±iT') = s(x(v)) (3.18) 

(here it is crucial that s is an homographic involution). For v in the strip ] — oo, oo[+i] — 
T',T'[ (which is mapped onto the complement of the cut [7_,7 + ]), let us define 

u{v) = x'(v)W hom (x(v)) . (3.19) 

Since Whom is meromorphic outside the cut [7_,7_|_], u is meromorphic on the strip. 
Note that ui is odd and 2T-periodic. Furthermore, by the above mapping property of 
x, we find that ui has well-defined values on the boundary of the strip, namely 

u(v±iT') = x'{v±\T') W hom (x(v±iT') =pi0), 0<v<T (3.20) 

(the value of u(v±iT') for other real values of v follows by oddness and 2T-periodicity). 
By (|3.13| ) and (|3.18| ) (and its derivative), we deduce that, for any real v, 

u(v + iT') + u(v - iT') - nu(v) = . (3.21) 

This relation implies that ui may be analytically continued into a meromorphic function 
defined on the whole complex plane, so that (|3.21|) holds for all complex values of v. 
For simplicity, we assume now n ^ ±2 (for n = —2, 2, the solution is also known, see 
for instance but we will not discuss it here). As discussed in [JT5| , any 2T-periodic 
meromorphic solution of ( |3.21|) can be expressed in terms of the "fundamental" solutions 
(±5, where rrb = arccos(n/2) and 

,,(01) M^Jk) iT 
S 1 ' 2T$! (-||r) (#|r) T 1 ' 

(see the Appendix in |T5[] for some properties of (b(v)). More precisely, any 2T-periodic 
solution of (|3.21|) is a linear combination of derivatives of translates of C±b- Here, we are 
looking for an odd solution with poles at ±Voo whose residues are fixed, when s(oo) ^ oo, 
by the condition 

4) = I>) ™ + ^M + .m ta„-«* (3.23) 



n 2 n + 2 x(v) J \v — v 

inherited from (|3.14 ), as well as poles at ± iT' with the same residues multiplied 
by n/2 = cos nb by virtue of ( ^.15| ). The general form for such a solution is 

u(v) = J2 a * {( {k) ( v ~ v °°) - - "oo)) (3-24) 

fc>0 
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where 

COO = Cb(V) \ C - b(V) , (3-25) 



and where the coefficients a& are determined by (|3.23| ). Note that the particular com- 



bination ( was chosen as it has residue 1 at and residue n/2 at ±iT' so that ( |3.23|) 
implies both ( |3.14j ) and ( |3.15| ). In the case s(oo) = oo (hence = iT'/2), the condition 
( P-23|) has to be modified according to (|3.16| ). We may observe that only a finite number 
of cifc's are non-zero when Vq is a polynomial. More precisely, when the maximal allowed 
face degree is D, then at = for k > D. 

After determining u in this way, it remains to satisfy the conditions p{p(±) = 0. By 
( [3.20Q , we may rewrite the spectral density ( 2.1 3|) in the parametric form 



/ / m /sx u(v + iT') — u(v — iT') m 

^ + lT » = 2i,4 + ir-) ■ < « < r ' < 3 ' 26 ' 

Since x has ramifications points of order 2 at iT' and T + iT' , the conditions p(7±) = 
amount to u(iT') = u(T + iT') = 0, which implicitly determine r y±. 

We further deduce that p (hence W) has generically a square-root singularity at 
7±, corresponding to a non-critical point. For T < oo, a critical point may only be 
obtained when uj" vanishes at either iT' or T + iT', so that p now decays generically 
with an exponent 3/2 at one endpoint of the cut: this corresponds to a generic critical 
point, in the universality class of pure gravity. Multicritical points can also be obtained, 
but only at the price of introducing negative weights in the model. 




— in 



Fig. 2: The (multivalued) change of variable from xtov, as given by ( |3.17|) 
in the particular case 7+ = s(7+), up some normalization factor which en- 
sures that v (7-) = ±i7r. The inverse function is given by Eq. (3.27). Here 
we assumed that s(oo) belongs to the cut ]7_, 7 + [ (solid line) so that the im- 
age (in the x-plane) of this cut by the involution s is ] — 00, s(7_)[U]7 + , 00 [ 
(dashed line). The cut maps onto the lines Imv = ±7r in the v- plane (solid 
lines), while its image maps onto the line Lm> = (dashed line). 



The only possible way to obtain a non- generic critical point with non- negative 
weights is to have 7+ = 5(7+), so that T = 00. The above strategy simplifies greatly in 
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this situation. Indeed, since one of the periods (2T) is infinite, the elliptic (doubly peri- 
odic) parametrization degenerates into a (simply periodic) trigonometric parametriza- 
tion (see Fig. 2). It is then convenient to rescale v by a constant factor so that 
i> (7—) = ±i7r (i.e. T' = 7r). The inverse function v 1— > x(v) then reads explicitly 

f \ 1 / w cosht> — 1 . . . „„ s 

«(») = (7 + -«(7-)) eaJlp _ eMh ^ +.(7-) (3-27) 

where is determined by the condition x(i7r) = 7_. Furthermore, for T = 00 and 
< n < 2, the fundamental solution ( in ( |3.25|) reads 

((v) = cosh(6f) cothf — sinh(6i>) (3.28) 

with 

irb = arccos ^— j . (3.29) 

Hence we may apply the above strategy on very explicit formulas. The only caution we 
must take concerns the condition p(7+) = which, since x'(v) ~ Ce~ v for v — > 00, now 
amounts to the condition u(v) = o(e~ v ). Noting that ( has the asymptotic expansion 

= e~ bv + e -( 2 - & )- + e -( 2 + & )- + {e-^ 2+b >) (3.30) 



with < b < 1/2, this yields a linear condition on the a^s appearing in fl3.24|) , that 
replaces the previous condition u(T + iT') = 0. 

Interestingly, this asymptotic expansion also yields the non-generic critical expo- 
nents. Indeed, we generically have oo(v) ~ C'e~^ 2 ~ b ^ v for v — > 00 and, noting that 
x{y) = 7-|- + ae~ v + o(e~ v ) with a > (since x(v ) approaches 7-1- from the right when u 
takes large real values), we deduce that Wh. m(x) ~ —a b ~ 2 C'(x — 7_|_) 1-6 . This is noth- 
ing but the leading singular term in the expansion of W around 7+, and the exponent 
1 — b is characteristic of the universality class of the dense 0(n) loop model on random 
maps. By transfer theorems, it amounts to the asymptotic behaviour 

i x r h C 7* 

r - -■ fc -)■ 00. (3.31) 



a 2 - h T(b-l) fc 2 - 6 ' 

Since < 6 < 1/2, the positivity of Fk (in a model with non-negative weights) requires 
that C > 0. When C = 0, we instead have u(y) ~ C"e-( 2+b)v so that W horn (x) ~ 
— a - ^ -2 ^^'^^— 7_|_) 1+& , and the exponent 1+6 is characteristic of the universality class of 
the dilute O(n) loop model on random maps. It corresponds to the asymptotic behaviour 

-ti +b c" 7 i 

Sr, ^^00 (3.32) 



a 2 + b r(-(6 + l)) fc 2 + b: 



hence C" < in a model with non-negative weights. It is not possible to have both 
C = and C" = in a such a model, thus ( p.31|) and ( p.32|) are the only possible 
"physical" non-generic critical behaviors. 
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4. The 0(n) loop model on cubic maps with bending energy: phase diagram 



4-1. The non-generic critical line 

In this section, we compute the non-generic critical line of the 0(n) loop model 
with bending energy defined in Section 3.2. As discussed in Section 3.3, non-generic 
criticality amounts to the condition 5(7+) = 7+. Note that 7+ is a continuous increasing 
function of the weights and h, with r y + = when the weights vanish. Then, when 
h is positive and a non-negative, this condition occurs when r y + reaches the smallest (in 
modulus) fixed point of the involution, namely 

Note that the other (largest in modulus) fixed point ^ a ^ h of the involution does not 
necessarily coincide with the other endpoint 7_ of the cut, but lies outside of both the 
cut and its image. The arguments of Section 6 shows that 7_ = ^ a \^ h can happen iff 
a = and the maps are bipartite. 

We now apply the strategy of Section 3.3 in the non-generic critical case. The 
trigonometric parametrization ( |3.27| ) reads here, for a 7^ 1, 

1 (1 — a) coshf + (1 + a) coshfoo 

x{v) = ^— , (4.2) 

(l — a z )h cosh f — cosh 

where 

„„ =aJgcosh Il^il±£|lzl. (4.3) 
(1 + a)(l + (1 - a)/i7_) 

The case a = 1 is obtained by taking a suitable limit, which we will discuss in the next 
section. Note that, since coshfoo is real, Voo necessarily lies on the boundary of the 
domain [0, oo[+i[0, tt]. 

Using this explicit parametrization, we may now determine the resolvent exactly, 
by computing the coefficients in ( |3.24| ). From now on, we will focus on the case where 
<7^ = 9$k,3, i-e. where faces not visited by the loops are themselves triangles, weighted 
by g. We have in this case Vq(x) = x — gx 2 so that ( |3.23| ) can be expanded into 

"{v) = Ta -^jz 1 -6 7- j u +2 



(4-n 2 )(l-a 2 )H {v- Voo y {v- Voo f 

{V-Voo) 2 (V-Voo) 
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with 



o 

8coth 
A* = 



A2 = 8 coth 2 Voo (a + sinh 2 Voo) 



B 2 = 4 coth 2 (1 - a 2 ) 

1 ' (4.5) 

A x = 4cothv °° ( 6 coth 4 Voo + (6a - 8) coth 2 Voo + 3(1 - a) 2 ) 
3 

Bi = 4 cothf QO (l — a 2 ) (coth 2 — (1 — a)) 
C = 2(l-a 2 ) 3 (n-2) . 



Then, since the fundamental solution £ in ( |3.28 ) has a simple pole at with residue 1 



satisfying (fO|) leads us immediately to the desired expression for a;: 

= (4-^1-^ 4 «> " "~> " C '"^ " "~» 

+ f + B 2 -^ ((> - ^00) - C"(-V - «„)) 

+ (Ax + (C'(v - Woo ) - C'(-« - «oo)) 

+ C (C(f - Uoo) - C(-« - u oo)) 



(4.6) 



In order for the solution (|4.6| ) to be consistent, we should ensure that the spectral density 
p vanishes at j±. As discussed in Section 3.3, demanding that p(j-) = amounts to 
the condition 

u(m) = (4.7) 



where, from 

cj(i7r) oc As-^ ^6cosh(6f oc ) cothvoo (b 2 + 6 sinh -2 v c 

- (2 sinh" 4 Voo + (36 2 + 4 coth 2 Voo) sinh" 2 Voo + b s ) sinh(fo 
+ + -^2-^ (6cosh(6woo) (2sinh" 2 Woo + 6) -coth^oo (6 2 + 2 sinh" 2 ^) sinh(ovoo) 



+ ( A i^ + S ipJ ^-sinh-^oo((6-2)sinh(^oo)-&sinh((6-2K 
+ Co sinh" 1 Vqo sinh((l — &)i>oo) • 

(4.8) 

Similarly, demanding that p(j+) = amounts, as already discussed, to the condition 

uj(v) = o(e~ v ) for v -)■ +00. (4.9) 
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By (|3.30|) and (^4.6|) , we have the asymptotic expansion 

u(v) = n(b)e- bv + k(2 - b)e- {2 ~ b)v + o{e- {2 - b)v ) (4.10) 



with 



«(6) 



(4 



_ n 2 ) | 1 _ a2)3 {- 4 3^(-26 3 sinh^ 0O )) + (a 2 ± + B2 l^j (26 2 cosh^ 00 )) 
+ (A^ + B^^j (-26sinh(6v oo )) + C (2cosh(?w oo ))} . 

(4.11) 

Then, the condition (|4.9|) reads nothing but 

«(6) = . (4.12) 
Finally, the positivity of the F^s requires by ( |3.31| ) the extra condition 

k(2-6)>0. (4.13) 

4-2. Phase diagram 

Fixing a > and b strictly between and 1/2 (i.e. < n < 2), the consistency 
conditions ([4.71) , (14-12 ) and (|4.13|) determine the values of g and h where the model 



is non-generic critical. Viewing Vqq cLS db parameter, (|4.7|) and (|4 . 1 2| ) form a system of 
inhomogeneous linear equations for the variables g/h 3 and 1/h 2 (or equivalently for 
the variables g/h and h 2 after multiplying the equations by h 2 ). Solving this system, 
we obtain a family of non-generic critical points, forming a line in the (g, h) plane 
parametrized by Voo- Recall that coshi^ is real, so that Voo varies on [0, oo[U(i[0, it]) U 
([0, oo[+i7r). This range is further restricted by the requirement that g > and h > 0, 
and by the condition (|4.13| ). We do not give here the precise parametrization of the 



non-generic critical line as the formulas are quite cumbersome, but we instead present 
in Fig. 3 the corresponding plots for b = 0.3 and several values of a. Plots for other 
values of b and a have essentially the same aspect: the line starts at g = and a finite 
value of h and ends at a point (g*, h*). 

All points on the line, except the endpoint (g*,h*), lie in the universality class 
of the dense 0(n) model: in particular has the asymptotic behaviour (|3.31| ) (with 
C = k(2 — b) > 0). The endpoint (g*, h*) lies instead in the universality class of the 
dilute 0(n) model: Fk has the asymptotic behaviour ( |3.32j ) (with C" = k(2 + b) < 0). 
On the plots, we have also indicated the existence of a line of generic critical points 
which links the point (g* , h*) to the generic critical point of pure triangulations without 
loops, obtained for h = and g = 1/(2 • 3 3 / 4 ). An explicit parametric representation 
for this line can be also obtained by the strategy of Section 2.5, but its expression is 



even more cumbersome than that of the non-generic critical line (see for instance |15| , 
Section 6.4] for its computation in the a = case, with different face weights of the 

form gf ] = g5 kA )- 
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Fig. 3: Phase diagram of the 0(n) loop model with bending energy in the 
(g, h) plane. The various plots correspond to b = 0.3 (n = 2cos7r&) and, 
from top to bottom, a = 0.5 (red), 1 (green) and 2 (blue). Plots for other 
values of b and a would essentially have the same aspect. The solid line 
corresponds to non-generic critical models and is computed exactly upon 
solving eqs.( f4.7|) and ( f4.12|) . In the absence of the third constraint (|4.13| ), 
this line would continue as indicated in dashed line (yellow). The condition 
( |4.13|) imposes that the line terminates at a point (g*, h*). All points on the 
non-generic critical line lie in the universality class of the dense 0{n) model 
except the endpoint (<?*, h*) which lies in the universality class of a dilute 
0(n) model. For the various models, a dot-dashed line was drawn (but 
not computed exactly here) to represent the line of generic critical points. 
This line connects the point (g*, h*) to the point ( 2 , 3 1 3 / 4 , 0) describing pure 
triangulations. The model is well-defined for value of g and h lying below 
the (generic or non-generic) critical lines. 



As already mentioned, the case a = 1 is special. Still, the equation of the non- 
generic critical line may in practice be obtained by simply taking the a — > 1 limit of the 
consistency conditions ( |4.7| ), ( |4.12j ) and (|4.13| ), with tending simultaneously to i7r/2 



as 



17T 

'>oo = TT + H a 



/i7_ 

2 ' 



+ 0((a-l) 2 ) 



(4.14) 



so that Eq. (fO|) remains satisfied. Solving the first two conditions in this limit leads 
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to the following parametrization of the non-generic critical line 



9_ 

h 

h 2 = 



2bV2Tnp-2V2 



n 



2bV2+7ip - v / 2^(l + p 2) 



bp 2 



6 V / 2T^(6 + (1 - b 2 )p 2 ) - 4(1 - 6 2 ) v / 2^np 



(4.15) 



48V4 - n 2 2bV2 + np- v / 2^ r n(l + ^P 2 ) 
where we use p = 1 — 2/ry_ = l- ^>0asa parameter along the line. The corre- 



sponding plot is displayed in Fig. 3. From the third condition, we find that the line 
ends at the dilute point 

y/(2-n)(6 + n) 



n 



(h 



*\2 



(2-b)b 4 



Q + n 12(l-6) 2 (2 + n)v/(2-n)(6 + n) ' 

These formulas are in agreement with the results of [fj| . In particular, we find the value 



h = i 

at this point. 



for the no n- generic critical point at g = 0, with 7_ 



(4.16) 
i value 



a =0.5 
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Fig. 4: Plots of the spectral density of the 0(n) loop model with bending 
rigidity at b = 0.3 and for a = 0.5, 1 and 2. For each case, we plotted the 
spectral density at the non-generic critical dilute point (g* , h*) (solid line), 
and at some (somewhat arbitrary) non-generic critical dense point along 
the non-generic critical line (dashed line). 
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For illustration, we have plotted in Fig. 4 the spectral density of eq. (|3.26| ) at the 
dilute point and at some (arbitrarily chosen) dense point on the non- generic critical line 
for b = 0.3 and for the three values of a displayed in Fig. 3. 

4-3. The n — > limit 




Fig. 5: Exact phase diagram of the 0(n) loop model with bending energy 
in the (g, h) plane, in the limit n — > and, from top to bottom, for a = 0.5 
(red), 1 (green) and 2 (blue). This phase diagram is similar to that of Fig. 3, 
apart from the fact that the line of non-generic critical points now has an 
infinite slope at its endpoint (g*, h*). We have g* = 2 ^ 3/4 independently 
of a, so that the generic critical points now form a vertical segment (dashed 
line). 

Fig. 5 shows the phase diagram obtained by our calculation in the n — > limit. In 
this limit, the equation for the non-generic critical line reads explicitly, in parametric 
form: 

g 2(a 2 - l)(coo + l)(oCoo +a - 1) 



h 



h 2(a-l) 2 +4a(a-l) Coo + (2a 2 + l)c4 

-4, 2(a-l) 2 + 4a(a-l)c 0O + (2a 2 -l)c 



(4.17) 



4(a 2 - l) 2 ( Coo + l) 2 2(a - l) 2 + 4a(a - l) Coo + (2a 2 + 1)4 



The parameter = coshfoo ranges a priori over the real line. By computing dis- 
criminants, we observe that the denominator in both fraction is always positive, while 
the numerator in the expression of h 2 has two reals roots at = — 2 ^^"^ . As a 

consequence, one should reject the values of Coo leading to negative values of f- or h 2 . 
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The line ends at the dilute point 



9 * 2-3 3 / 4 ' k * (1 + a) 3 3 / 4 + 3V4 (4 ' 18) 

and, since the value of g* matches that of the (generic) critical point of pure tri- 
angulations, we expect that the generic critical line becomes the vertical segment 
{{g*,h),0<h<h*}. 

Let us now explain how to recover this phase diagram without recourse to our 
calculation by the following simple argument. The n — > limit corresponds to config- 
urations of maps with a fixed finite number of loops. The generating function for such 
configurations may be expressed in terms of the generating function for pure triangula- 
tions with a boundary, gathered in the resolvent. For pure triangulations, this resolvent 
has a cut with endpoints S ± 2\/R where R and S satisfy the relations (see for instance 
0) 

R = l + 2gRS, S = g(2R + S 2 ). (4.19) 

Now criticality may be achieved in two different ways. The loops may remain 
small, in which case a generic critical point is simply obtained whenever dR/dg = oo, 
i.e. g = g* = 1/(2. 3 3,/4 ) irrespectively of h (small enough). We therefore deduce the 
existence of a vertical line of generic critical points in the (g,h) plane at g = g* . On 
the other hand, non-generic criticality occurs when the loops become large. Now a loop 
with contours of lengths k and k' yields a contribution of order F£ ure tn pP^ Ye tri Ak y k' to 
the partition function. In particular, loops of length K (i.e. with k + k' = K) contribute 
at large K as 

(S + 2^fR) K A W -{S + 2^R) K ((a + l)h)f (4.20) 

k+k'=K 

(the reader may easily derive the identity ^2j~+k'=K A k,k' = i x + +xZ K )/2 where x± = 
l/((a± l)h) are the fixed points of the involution). The condition for the occurrence of 
large loops reproduces precisely the non-generic critical condition ( |4 . 1| ) , namely 



7+ 



(a + l)h 



S + 2VR. (4.21) 



Using ( f4.19|) , this leads immediately to the following parametric form of the non-generic 
critical line: 

| = (1 + a) (VMl " o) + °) 

2 = a(l-a)(l-2a) (4 ' 22) 
9 2 

where we have set a = gS (this parametrization has been introduced in |T£|] ). The value 
of the parameter o varies along the line between and a* = (3 — v3)/6, corresponding 



to the dilute point (|4.18[) , which is also the endpoint of the vertical line of generic critical 
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points. Now it is a straightforward exercise to check that the parametrizations ( [4.17 ) 
and (|4.22|) are indeed fully equivalent upon relating the varying parameters a and 
through 

= 2(q Coo + a- l) 2 

2(a - l) 2 + 4a(a - l) Coo + (2a 2 + 1)4, " 1 " ' 



5. The twisting loop model on quartic maps 

The analysis of Section 4 can be applied to obtain the non-generic critical line of the 
twisting loop model defined in Section 2. To be explicit, we decide to set = gSj-^, 
i.e. the faces not visited by the loops are themselves squares, weighted by g. Eq. ( |2.26| ) 
reduces in this case to 

W(X + iO) + W(X -iO) = l-gX-2n W^ 1 - X) , (5.1) 

of which we look for a solution W(X) having a cut along the real segment [0, T] for some 
positive T (recall that the resolvent W(x) is recovered via W(x) = x W(x 2 ) and has a cut 
along the real segment [7_, 7+], with 7± = ±vT). Equation (|5.1| ) is similar to (yet even 
simpler than) ( |3.7| ) and we may solve it by the strategy of Section 3.3. It now involves the 
particular involution X 1— > s(X) = h^ 1 — X, here n is replaced by 2n and the polynomial 
Vq' reads Vq(X) = l-gX. As before, we write W(X) = W part (X) + W hom (X), with 
the particular solution (with no cut) 

W (X) ( 1 -9X)-n(l-g(h^-X)) 

W^ t {X) - 2(1 -n 2 ) (5 ] 

and where W^o^X) is now a solution of the homogeneous equation 

W hom (X + iO) + W hom (X - iO) + 2n W hom (h^ - X) = . (5.3) 

Again, among all solutions of this equation, W / hom(^) is determined by requiring that 
W(X) has no poles at finite values of X and that W(X) ~ 1/X for X — > 00. We are 
here in a situation where s(oo) = 00 so that, as for loops on triangulations with a = 1, 
we have a single condition to satisfy, namely 

W, (X)- £ ~ 9X) - n(l - gfci - X)) 1 / 1\ 

w hom (x)- ^r^) + x +0 {x^) ■ (5 - 4) 



The general solution of ( |5.3| ) is still given by elliptic functions. Note however that 
in the elliptic integral (|3.17|) , 7+ and 7_ should be respectively replaced by V and 



0. We again concentrate on the non-generic critical case which occurs whenever the 
endpoint V of the cut is a fixed point of the involution s(X), namely V = 1/(2/12). The 
parametrization Q3.27|) is here replaced by 

X(v) = ^-(l + -^—]. (5.5) 



2h<2 V coshf 
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Defining 

u(v)=X'(y)W bom (X(v)) (5.6) 
as in ( p.!9|) , we find that u(v) is an odd meromorphic function satisfying the functionai 



equation 

u>(v + br) + u(v — br) — 2nu(v) = , (5-7) 



while Eq. ( |5.4j) translates into 

9 1 i i(jg-2h 2 ) I I 



UJIV 



8/jl(l-n)( v -f)3 8fel(l + n)( v -f) 



2 

Again, cD(t>) may be expressed in terms of the fundamental solution ( of Eq. ( |3.28|) , but 
now with b given by 

%b = arccos n (5.9) 
which lies in the range < b < 1 /2 for < n < 1 . We get explicitly 



8fel(l-n 2 ) V V 2 ! s V 2 

1 'dv-%\-i' ' m 



l + Ti v V 2 / v 2 

Since W(x) = x W(x 2 ), the spectral density is recovered via 

p(x) = M (v^ hom (x 2 + iO) - V^ hom (x 2 - iO)) (5.11) 

so that 

P (±V^Tto)) = ^2±M . "(» ± ^) ; gfe - !^ (5 . 12 ) 

V / + 17T) 2l7T 

for real positive v. Contrary to the previous case, we need not have to(m) = 0, since 
p(0) does not vanish. But we must impose again the condition u(v) = o(e~ v ) for v — > oo 
so that p{y±) = 0. By ( |3.3U| ) and ( |5.10| ), we have the expansion 

u(v) = k(b)e- bv + k{2 - b)e-( 2 ~ b > + o{e- {2 - b)v ) (5.13) 



with 



_ (g 6 2 + 16 ft| (n - 1)) cos f - 2 6 (g - 2 fe 2 ) sin f 
KW " Sh 2 (l-n 2 ) • (5 - MJ 



The condition k(b) = yields 



&= 4^-16V^M. (5.15) 
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Finally, as before, the positivity of requires that h(2 — 6) > 0, which yields the 
condition 



9< 



4(2-6) h 2 + 16V1 - n 2 h\ 



(2-6) 2 + (2-6), /t 



l+n 



The non-generic critical line therefore ends at the dilute point 

6(2 - 6) 



9 



4(vT^I+ (1 - ^v/TT^) 2 

6(2 - 6) 

8(1 -n+ (1 -6) VI -n 2 ) ' 



(5.16) 



(5.17) 



/7 
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Fig. 6: The phase diagram of the twisting loop model in the variables 
(g,h2) for 6 = 0.3 (n = cos7r6). The line of (dense) non-generic critical 
points (solid red line) is the portion of the parabola (|5.15|) (dashed yellow) 
delimited by the conditions g > and ( J5.16| ). The line ends at the dilute 
point (g*,^) of eq. ( |5.17| ). We indicated in dot-dashed line the existence 



of a line of generic critical points linking the point (g* 
^1/12,0) describing pure quadrangulations. 



h 2 ) to the point 



The phase diagram of the twisting loop model for 6 = 0.3 is displayed in Fig. 6. 
A line of generic critical points links the dilute point (g*,/^) to the critical point of 
pure quadrangulations without loops at g = 1/12, h% = 0. For illustration, we have 
also plotted in Fig. 7 the spectral density at the dilute point and at some (arbitrarily 
chosen) dense point on the non-generic critical line. 
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dilute 
dense 



Fig. 7: Plots of the spectral density of the twisting loop model at b = 0.3. 
We plotted the spectral density at the non-generic critical dilute point 
(g*, h%) of eq. ( p.l7| ) (solid line), and at some (arbitrary) non-generic crit- 
ical dense point along the non-generic critical line (dashed line). 




1/12 



Fig. 8: Phase diagram of the twisting loop model at n = 0. The line of 
non-generic critical points now has an infinite slope at its endpoint (g* , h^). 
We have g* = 1/12 so that the generic critical points now form a vertical 
segment (dashed line). 



When n — > 0, the dilute point (p.!7| ) has coordinates g* = 1/12 and h\ = 1/16, with 
the same value of g as for critical pure quadrangulations, so that the generic critical 
line becomes a vertical segment {(1/12, /i), < h < 1/16} (see Fig. 8). As for the 
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equation for the the non-generic critical line, it may be recovered by simply writing 
that T = (7+) 2 is the fixed point l/(2h 2 ) of the involution, with, for n — > 0, a value of 
7_i_ equal to that of pure quadrangulations. More precisely, we have as in Section 3.5 
7+ = S + 2v^R, but now with S = and R solution of (see for instance |f[8| ) 

R=l + 3gR 2 . (5.18) 

Writing (2y/R) 2 = l/(2h 2 ) leads to 

9=l(h 2 -8h 2 2 ) (5.19) 
which matches ( |5.15| ) specialized to n = 0, i.e. 6 = 1/2. 



6. Analytic properties of the resolvent for maps with controlled face degrees 

In this Section, we review a number of analytic properties of the resolvent for maps 
with controlled face degrees, as we have used them in Section 2 in the framework of our 
loop models. 

6.1. The one-cut lemma for maps with controlled face degrees 

We say that a sequence of non-negative numbers (gk)k>i is admissible when the 
generating series of pointed rooted planar maps with weight gk per face of degree k is 
finite. In particular, a model with bounded face degrees corresponds to having all but 
a finite number of gkS non-zero, and such a sequence is admissible when these weights 
are small enough. 

We denote as before by Tpigi, g 2 , ■ ■ •) the multivariate generating function for maps 
with a boundary of length p (with the convention J-'q = 1), which is finite whenever the 
sequence (gk)k>i is admissible, and we introduce the resolvent 

Wfr) = (6.1) 

p>0 

(the term resolvent, and the convention of taking W(x) as a series in 1/x are borrowed 
from random matrix theory). The following property is essential to the determination 
of W(x) in specific models: 

• The one-cut lemma: When (gk)k>i is an admissible sequence, there exists a real 
segment / = [7_,7_|_] such that x i— > W(x) defines an analytic function on C \ /, which 

has a finite discontinuity on /. Moreover, we have | — y | < 7+ with equality iff g\ = for 

all odd I. Eventually, the spectral density p, defined for x E I by 

„(*) = w(»-io)- >v(» + io) , (6 . 2) 
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vanishes at the endpoints of /, and is positive in the interior of I. 



This lemma essentially follows from the discussion in |H| for the case of a model with 
bounded face degrees. In this paragraph, we extend it to the case of non-necessarily 
bounded degrees by describing the discontinuity set of W. We use the relations estab- 
lished in [^0[ thanks to the bijection between maps and mobiles, and the relation of the 
latter to Motzkin paths. We need some notations: first, let J- p [u] (resp. ^'[w]) be the 
generating series of rooted (resp. pointed rooted) maps with root degree p where, on 
top of the face weights, each vertex of the map is counted with a weight u. Since the 
two definitions differ only by the marking of a vertex, we have 

<LF p [u] , , 



In other words 



-1 du 

F P (gi,g 2 ,...)= F P [u=l]= I —J*[u] . (6.4) 



o 



Then, let P P (R, S) be the generating series of lattice paths in Z 2 from (0,0) to (p,0), 
where the set of allowed steps consists of: 

- level-steps — > (i + counted with weight S; 

- up-steps — > (i + l,j + 1) counted with weight a/R; 

- down-steps — > (i + 1, j — 1) also counted with weight \/R. 
An easy counting argument shows that 



LP/2J 



w*)=Z(^*- s - («■<» 

We may pack those numbers in the generating series 

p( x R C) = P p( R > S ) = 1 (Q Q) 

^ xP+1 Vx 2 -2Sx + S 2 -4R ' V ' ; 



We observe that there exists a segment I(R, S) = [S — 2yR, S + 2yR] such that x 
P(x, R, S) is analytic in C \ I(R, S) and has a discontinuity on I(R, S). 

The relation established in |20|1 is [u] = u P p (R[u], S[u]), where R[u] and 

S[u] are the unique power series in u and the gkS determined by the system 



S[u]=Y,g k P k - 1 (R[u],S[u]) 
R[u] =u-^p- + ^J29k Pk(R[u],S[u]) 



2 

fe>l 



Actually, R[u] and S[u] have themselves combinatorial interpretations in terms of maps, 
and 2R[u] + S 2 [u] is the generating function for pointed rooted maps. In particular, 
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R[u] and S[u] are continuous increasing functions of u G [0, 1] whenever the sequence 
(9k)k>i is admissible, with R[0] = S[0] = 0. From (|6.1|) and (|6.4|), the resolvent can be 
written as 

W ^ = f 7=2 ocr i tc 2 r i m i ■ ( 6 - 8 ^ 
Jo V x 2 - 2S[u] x + S 2 [u] - AR[u] 

We deduce that W(x) is an analytic function outside the discontinuity set 

1= |J [5M-2^M, i S'M+2/RM]. (6.9) 
ue[o,i] 



It is clear that u S^ia] + 2yR[u] is a strictly increasing function of u whose derivative 
does not vanish. Although it is not obvious, we justify in Section 6.3 that u h->- S[u] — 
2vAR[m] must be strictly decreasing with a non- vanishing derivative. Hence 

I=[ T -,7+] (6-10) 

where 

7 _ = S[l] - 2y/W], y + = S[l]+2y/R^] (6.11) 

In particular, we have |-y | < 7+, with equality if and only if S[l] = 0. From (|6.7|) , this 

is equivalent to having g\ — for all odd integers /, which means that only faces of even 
degree are allowed. By (|6.8|) and by monotonicity of u 5 [it] ± 2^/i?[u], the spectral 
density reads, for x g]7_,7+[, 

= ~ f 7m i cart ocr i ^ 
K Ju(x) ^AR[u] - S 2 [u] + 2S[u] x-x 2 

where u(x) is obtained via S[u(x)] + 2yJ R[u(x)] = x for x > and via 5'[w(x)] — 
2yJ R[u{x)\ = x for x < 0. Clearly, p is positive on ]7_,7_|_[ and the non- vanishing 
of (d/du)(S[u\ ± 2^jR[u\) ensures that p is continuous and vanishes as a; — > 7±, as 
announced. 



6.2. Functional relation for maps with controlled face degree 

We now assume that (gk)k>i is an admissible sequence, so that the one-cut lemma 
applies. The discontinuity of W on its cut is related to the spectral density p by ( |6.2|) 
and conversely, W is the Stieltjes transform of p, namely 

= P + p^y)dy^ = r+ yk 

We claim that the resolvent W(£) is characterized by W(x) ~ 1/x and by 

VxG]7_,7+[ W{x + iO) + W(x - iO) = V\x) , (6.14) 
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where V is the so-called "potential" defined via 

V\x) =x-Y,9kX k ~ 1 . (6.15) 

fc>i 

In terms of the spectral density, this characterization is equivalent to Jj p(x)dx = 1 and 

p(y)dy V'(x) 



VxG]7_,7+[ p - v -/ 



x-y 



(6.16) 



In this paragraph, we shall justify this functional relation for general admissible se- 
quences. 



Tutte's recursive decomposition |2T| of rooted planar maps allows to write, for any 
sequence of face weights, 



m-2 



Fm = ^ Fk?m-2-k + '^29kJ r m+k-2 ■ (6-17) 

k=0 k>l 

If we want to write this relation in terms of the resolvent W(x), we first need to know the 
location of the singularities of V'(x) which, at first, only makes sense as a power series 
in x. From (|6.7|) and the fact that S[l] and R[l] are finite for admissible sequences, we 
can deduce a lower bound on the radius of convergence of V'(x). We use the asymptotic 
behavior of S[l]) when k — > oo, which can be extracted from ( |6.6| ) by transfer 

theorems. Since |7_| < 7+, 7+ is always the dominant singularity of P(x, R[l], S[l]) 
and P k (R[l], S[l]) = 0(k~ 1/2 'j^ k ) (when S[l] = 0, 7 _ = -7+ and P k vanishes for odd 
k but the same estimate holds for even k). Accordingly, the radius of convergence of V 
is at least 7 + . In particular, x h-> V'(x) defines a real-analytic function at least on the 
segment ]7_,7+[. So, we may represent, for any integer m > 1 

X T" T* "D /-Tfl—1 ! z \ /-k—1 ^~ m-\-k — 2 

— Ok,2-Tm + 9k J-m+k-2 = KeS^ \Ok,2 ~ 9k)Q 



£ (m+fc _ 2 ) + l 



Resr- L (Sk,2-9k)e- L mo re.is) 



™ Jc 



where the contour C surrounds the cut [7- , 7+] in positive orientation, and is included 
in the domain of analyticity of V . To be precise, if V happens to have a singularity 
at £ = 7 + , we choose a contour sticking to this point when going around the cut. The 
same caution might be necessary at £ = 7_ = —7+ in the case of maps with even face 
degrees. This allows to do the summation over k, namely 

V -5k, 2 J C m + gkJ C m+k-2 = -7^- <f <% ^ ^'(OMO . (6.19) 

kTi 2l7r J c 
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Now, let us take x G C such that \x\ > 7+, and sum over non-negative integers m 
the equations (|6.17| ) with a weight x~ m . All the series involved are convergent and we 
obtain 



W 2 (x) - J- / V'(£) W(£) = . (6.20) 

This is an equality between functions of x which are analytic in the unbounded connected 
component of C \ C, so it must be valid for any x in this domain. In particular, we can 
squeeze C to stick to the cut /, and compute the discontinuity of the left hand side of 
( |6.20| ) at a point x interior to /, which should equal zero. We find the desired result 

VxG]7_,7+[ W{x + iO) + W(x - iO) - V'{x) = . (6.21) 

When |7_| < |7 + |, this relation also holds at x = 7_. 

6. 3. Edge behavior of the spectral density 

We now justify that u h-> S[u] - 2y/R[ u] is a decreasing function of u with a non- 
vanishing derivative. We may eliminate the case S[u] = (i.e. we assume that at 
least one gi for odd I is non-zero) as the property is obvious in this case. Our proof, 
as presented here, relies on real analysis rather than combinatorics. We introduce the 
shortcut notations T[u] = ^/R[u], which is a strictly increasing function of u, and 
^(0 = 12k>i 9k By the same trick that we used in Section 6.2, the relations 

can be rewritten 



2i7r J c ^/|3 - 2S[u}£ + S 2 [u] - 4T 2 [u] 
2T [u\ + S [u\ — 2u-\ 

2l7T J c 



(6.22) 



T c y /^_2S[u]C + S 2 [u] -4T 2 [u] 



These relations are also well-known, and we have justified them for arbitrary (non- 
negative) admissible weights. 

If we use the change of variable £ = S[u] + 2T[u] cosy? in the integrals, we find 



S[u] = - I dip U (S + 2T[u] cos <p) , 
n Jo 

T 2 [u]=u+[ — / dip cos (p U(S[u] + 2T[u] cos cp) )T[u] . 
V7r Jo ' 



(6.23) 



We then differentiate these relations with respect to u, and use an integration by parts 
to arrive at 

S\u] = S[u] L Q [u] + 2T[u] Li[u] , 

(6.24) 

2T[u]T[u] = 1 + 2T[u]T[u] L [u] + T[u]S[u] L^u] , 
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where we have set 



L Q [u] = - [ dipU'(S[u} + 2T[u] cos^) , 

7T Jo 

i r 

L\[u\ = — dtp cos ipU (S[u] + 2T[u]cos(p) 
Jo 



(6.25) 



Eqs. (|6.24[ ) could have been obtained alternatively upon differentiating directly ( |6.7| ) 
with respect to u. The quantities Lq[u] and Li[u] would then be identified as generating 
functions for paths with height difference or 1 at their endpoints, and with a marked 
point along the path. As such, Lq [u] and L\ [u] are clearly positive. In our setting, this 
property follows simply from the expressions 



(6.26) 



with T[u] > 0, while all derivatives of £ i— > U(£) at S[u] > are positive since all 
the gk's are non- negative (here we eliminate the trivial case where only gi would be 
non- vanishing) . Combining the relations ( |6.24| ) and eliminating Lq[u], we deduce 

™) 2 "(%]) 2 = ^L>0. (6.27) 



Since S[u] and T[u] = y/R[u] are strictly increasing functions of u, the claim follows: 
S[u] - 2T[u] < 0. 

6.4- Remark on complex valued face weights 

Let us say a word on more general, maybe complex valued face weights. To avoid 
confusion in this paragraph, we shall denote by 1Z and S the power series in g^ satisfying 
(|6.7|) . We say that a sequence of complex numbers (gk)k>i is sub- admissible if the 
sequence (|<7fc|)fc>i is admissible. In this case, lZ[u] and S[u] assume finite values for the 
sequence (|<7fc|fc) and any u G [0, 1], that we denote Rm[u] and Sm[u}- Hence, they also 
assume finite values that we denote R[u] and S[u] for the sequence (gk)k>i itself, and we 
have \R[u}\ < Rm[u] and \S[u]\ < Sm[ u ] f° r all u. Moreover, the radius of convergence 
of V is at least 7+,m = ^Affl] + 2a/-Rm[1]- Recalling that Sm[w] + 2^Rm[u] increases 
with u, we have \S[u] ± 2y/ R[vi\ \ < 7+,m- Thus, the disk of convergence of V contains 
the interior of the discontinuity set of W(x) computed with the weights (gk)k>i- This 
ensures also for sub-admissible sequences the validity of ( |6.20|) . The discontinuity set is 
now: 

/= |J{ xG C (x 2 - 2S[u]x + S 2 [u] -AR[u]) < 0} (6.28) 
we [o,i] 
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It can be quite complicated, since R[u] and S[u] moves in the complex plane when u 
increases. Nevertheless, we find: 

W(x+) + W(x-) -V'{x) = (6.29) 

for any point x of / such that \x\ < 7+,m arid such that the intersection of / with a small 
ball B centered at x is a piecewise C 1 curve separating B in two connected components 
B + and B_ which both intersect the unbounded component of C \ /. x then denote 
a sequence of points in B± which converges to x. 



7. Conclusion and discussion 

In this paper, we have shown how the gasket decomposition, upon relating loop 
models to models of maps with controlled face degrees, immediately provides a func- 
tional equation for the resolvent of the model. We then looked for the largest possible 
class of models for which this functional equation may be solved by the same techniques 
that were used to solve earlier models. In this respect, it was crucial that the cut py_, 7+] 
interacts with a single mirror image, which is guaranteed if the ring (grand canonical) 
generating function has a single pole, in which case a homographic involution naturally 
enters the game (if the model is symmetric). For non- negative weights, this involution 
is moreover locally decreasing but the solving technique could easily be extended to deal 
with increasing homographic involutions. We discuss in Appendix A how a non-generic 
critical point could emerge in this case. 

All the above solvable models (with non-negative weights) can be realized within 
the framework of a single model of loops with bending energy, or for twofold avatars 
of such models in bipartite maps (such as the twisting loop model). All those models 
display the same critical behaviors as those observed in previously solved 0(n) models 
and no new universality class is found. In particular, at non-generic critical points, the 
critical behaviors are also those of maps with large faces, namely the stable maps of 
II- 

The problem of finding a resolvent satisfying (|2.22|) for the quadrangular case with 



arbitrary hi and /12 is much more involved as the cut now interacts with two mirror cuts 
(its images by y+ and y_). This problem is still open and may lead to new interesting 
critical phenomena. In all generality, loop models may involve several poles for the ring 
generating function and new techniques are needed to solve these multi-pole models. 

As a more promising issue, it is worth mentioning that the nested loop approach 
can be generalized to describe Potts model (possibly with defects). This again leads to a 
functional equation and it may easily be seen that it now involves a pair of homographies 
which are inverse of one another, rather than a single involution. This problem is under 
investigation. 
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Appendix A. Case of a general nomographic involution 

We have seen that a model with non-negative local weights cannot lead to a homo- 
graphic involution s(x) which is increasing, or for which s(0) < 0. Let us nevertheless 
consider here the case of an arbitrary homographic involution, namely 

o_-|^ a)jMGR . (A.l) 

p — ox 

Recall that s has real fixed points iff it is locally decreasing over R, i.e. iff (3 2 — ad > 0. 
Such a general homographic involution may be obtained in loop models for a particular 
set of real weights. We shall assume here that these weights may be chosen in such 
a way that the resolvent has one cut only on a segment I = [7- , 7+] of the real line. 
This requirement is not an empty condition since the one-cut property is not automatic 
when the weights are not all positive. Note however that, under this assumption (and 
provided the weights are real), the spectral density must be positive in the interior of 
I. 

We may distinguish six situations, displayed in Fig. 9, depending on the local mono- 
tonicity of s and on the relative position of I and s(I). We shall now examine these 
situations and explain when a non-generic critical behavior can be reached. 

The situation 1— corresponds to a locally decreasing involution with s(I) > I. The 
non-generic critical behavior is obtained when 7 + collides with 5(7+), i.e. is one of the 
fixed points of s. This case can always be realized in a model with bending rigidity (see 
Section 3.2), and we have already studied it in details. 

The situation 2— corresponds to a locally decreasing involution for which s(J) < I. 
The non-generic critical behavior is obtained when 7_ collides with s(7_), i.e. is again 
one of the fixed points of s. 

The situation 3— corresponds again to a locally decreasing involution for which 
the pole of s belongs to I so that s(I) is split into two semi-infinite intervals. The non- 
generic critical behavior is obtained when 7+ collides with 5(7+), or when 7_ collides 
with s(7_), or both. This case also arises in a model with bending rigidity. 

The situation 1+ (resp. 2+) corresponds to a locally increasing involution for 
which s(J) > I (resp. s(J) < I) and does not feature any non-generic critical behavior. 
Indeed, if r y± collides with s(7=p), by involutivity both extremities of I should collide 
simultaneously with the extremities of s(J), which is clearly not possible. 

The situation 3+ corresponds to again to a locally increasing involution for which 
the pole of s belongs to I, so that s(J) is split into two semi-infinite intervals. A non- 
generic behavior is obtained whenever 7 + collides with s(7_), which implies that 7_ 
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Fig. 9: The six situations for the relative positions of / and s(J) (see text) 



collides with 5(7+) by involutivity, so that both extremities of I collide simultaneously 
with the extremities of s(I). 

The techniques presented in this article to determine the resolvent W(x) can easily 
be applied to the six situations. The equations of the critical manifold, and in particular 
the expression of W(x) on the non-generic critical manifold in terms of trigonometric 
functions would follow. 
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